Weak and strong measurement of a qubit using a switching-based detector 
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We analyze the operation of a switching-based detector that probes a qubit's observable that does 
not commute with the qubit's Hamiltonian, leading to a nontrivial interplay between the measure- 
ment and free-qubit dynamics. In order to obtain analytic results and develop intuitive understand- 
ing of the different possible regimes of operation, we use a theoretical model where the detector is a 
quantum two-level system that is constantly monitored by a macroscopic system. We analyze how 
to interpret the outcome of the measurement and how the state of the qubit evolves while it is being 
measured. We find that the answers to the above questions depend on the relation between the 
different parameters in the problem. In addition to the traditional strong-measurement regime, we 
identify a number of regimes associated with weak qubit-detector coupling. An incoherent detector 
whose switching time is measurable with high accuracy can provide high-fidelity information, but 
the measurement basis is determined only upon switching of the detector. An incoherent detector 
whose switching time can be known only with low accuracy provides a measurement in the qubit's 
energy eigenbasis with reduced measurement fidelity. A coherent detector measures the qubit in its 
energy eigenbasis and, under certain conditions, can provide high-fidelity information. 



I. INTRODUCTION 

The quantum theory of measurement is of particular 
significance in the study of physics because it lies at the 
boundary between the quantum and classical worlds 
Quantum measurement has received renewed attention 
recently from a new perspective, namely its role in quan- 
tum information processing 0. Understanding how the 
quantum state evolves during the measurement process 
and the information that can be extracted from the raw 
measurement data is crucial for optimizing the design 
and operation of measurement devices. 

One typically thinks of the measurement process in 
terms of the instantaneous, projective measurement en- 
countered in the study of basic quantum phenomena. 
For example, one sometimes imagines that at any de- 
sired point in time the spin of a particle can be measured 
along any desired axis. In many situations, however, the 
measurement device is weakly coupled to the quantum 
system being measured, such that the latter evolves al- 
most freely according to its free Hamiltonian while being 
barely perturbed by the measuring device. 

In the context of solid-state quantum information pro- 
cessing, the weak coupling between the qubit and the 
measuring device is imposed by the need to minimize 
the decoherence of the qubit before the measurement; 
in practice it is not possible to completely decouple the 
qubit from the measuring device, even before the mea- 
surement is performed. It is worth mentioning here that 
although this constraint of weak coupling might seem to 
be undesirable for measurement purposes, weak coupling 
between the qubit and the measuring device can have 
advantages over strong coupling. For example, when the 
energy eigenstates of a multi-qubit system contain some 
residual entanglement between the qubits, this entangle- 



ment translates into errors when using strong single-qubit 
measurements, whereas weak coupling to the measuring 
device can lead to a measurement that is protected from 
such errors [3]. 

In the study and implementation of qubit-state read- 
out, one strategy to remove any conflict between the 
qubit's free dynamics and the measurement process has 
been to bias the qubit at a point where the qubit Hamil- 
tonian and the operator being probed by the measuring 
device commute. In this case the free qubit dynamics 
conserves the quantity being measured. This dynamics 
can therefore be ignored during the measurement pro- 
cess, and the only difference between weak and strong 
qubit-detector coupling is in the amount of time needed 
to obtain a high-fidelity measurement result. It should be 
noted, however, that in practice it is not always possible 
to reach such an ideal bias point. 

Here we are mainly interested in the analysis of the 
measurement process when the qubit Hamiltonian and 
the operator being probed by the measuring device do 
not commute. One must therefore analyze the simulta- 
neous working of two physical processes that, in some 
sense, push the qubit in different directions. We shall fo- 
cus on the experimentally relevant example of a binary- 
outcome switching-based measuring device (i.e. the sit- 
uation where the measuring device can switch from its 
initial state, to which we refer as 0, to another state, 
called 1, depending on the state of the qubit), mainly as 
applied to superconducting circuits Our work has a 
similar motivation to the recent experimental and theo- 
retical work of Refs. [1, 0]. However, that work consid- 
ered a measurement protocol where the measured quan- 
tity is the value of the bias current at which switching 
occurs, whereas we consider the more common approach 
that uses the switching probability at fixed bias param- 
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eters. Other related recent studies include work on the 
continuous monitoring of Rabi oscillations 0] , as well as 
the theoretical study in Ref. Q of multiple linear detec- 
tors simultaneously probing different observables of the 
same qubit. A models similar to ours was also used in 
Ref. [9| to study certain aspects of quantum state evolu- 
tion under the influence of weak measurement. 

It is worth clarifying from the outset a point concerning 
terminology. A strong measurement is a measurement 
whose outcome is strongly correlated with the measured 
quantity, thus specifying with a high degree of certainty 
the state of the measured system. After the measure- 
ment the system is, to a good approximation, projected 
onto a state with a well-defined value of the measured 
quantity. A weak measurement, on the other hand, is a 
measurement whose outcome is weakly correlated with 
the measured quantity, thus serving as a weak indicator 
of the state of the measured system. Under suitable con- 
ditions a sequence composed of a large number of weak 
measurements can result in an effectively strong measure- 
ment. A weak measurement can also result in reliable 
information when repeated a large number of times on 
identically prepared systems. It is sometimes stated that 
weak measurement results from weak coupling between 
the measured system and the measuring device. It should 
be noted, however, that weak coupling does not necessar- 
ily imply a weak measurement. We shall therefore use the 
terms slow-switching and fast-switching detectors in this 
paper. This distinction could be understood as being be- 
tween weak-coupling and strong-coupling measurements. 
A weak-coupling measurement can result in a weak or 
strong measurement depending on other system parame- 
ters, as we shall show with a number of examples below. 

This paper is organized as follows: In Sec. II we de- 
scribe the model and briefly explain an example demon- 
strating the operation of the switching-based detector in 
a simple case. In Sec. Ill we outline our approach to the 
analysis of the measurement process when using a detec- 
tor that has a short coherence time. A number of specific 
examples are treated in Sec. IV. The case of a detector 
that has a relatively long coherence time is treated in 
Sec. V. The results of Sees. II-V are reviewed in Sec. VI. 
We discuss the relevance of our analysis and results to 
the readout of superconducting flux qubits in Sec. VII, 
and we summarize our results in Sec. VIII. 



II. MODEL 

We consider the situation where one is trying to read 
out the state of a quantum two-level system, i.e. a qubit. 
The qubit is placed in contact with a binary-outcome 
detector, i.e. a detector that produces one of two possi- 
ble readings. The detector is initialized in one of its two 
possible outcomes, and it can switch to the other one, 
depending on the state of the qubit. For a large part of 
the following analysis, we shall assume that the detector 
itself is a quantum two-state system that is constantly 
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FIG. 1: Schematic diagram of the theoretical model that we 
use to analyze the operation of a switching-based detector. 
The qubit's energy eigenstates are denoted by |0) and and 
they have energy separation E. The detector is a two-state 
system (TSS) [with a non-switched state lA'^) and a switched 
state 15*)] that is constantly being monitored by a macroscopic 
system. The detector probes the operator of the qubit. x 
represents the detector's operator through which the detector 
is coupled to the qubit. The macroscopic system monitors 
whether the TSS remains in the initially prepared state |A'^) 
or switches to the state IS). 



being monitored by a macroscopic, classical system (see 
Fig. 1). We shall refer to the non-switched and switched 
states of the detector as |A^) and |S'), respectively. Al- 
though the two-state-detector model does not accurately 
describe the realistic measuring devices that we shall dis- 
cuss in Sec. VII, this simplification in the model will allow 
us to carry out our calculations analytically and obtain 
a good amount of insight into the dynamics of the mea- 
surement process. We shall come back to this point and 
discuss the limitations of the simplified model in Sec. VII. 

The detector is initially prepared in the state \N). 
We shall assume that the switching is irreversible, such 
that after a switching event the detector remains in 
the switched state until it is reset by the experimenter. 
One possible way to naturally prevent the detector from 
switching back to the state \N) is to replace the state \S) 
by a continuum of states. However, this change would 
results in lengthier algebra and similar results to the 
two-state-detector model. We therefore do not make this 
change, and we simply impose by hand the requirement 
that once the detector is observed to be in the state \S) 
it does not switch back to the state \N). For most of this 
paper, we shall also ignore any intrinsic qubit decoher- 
ence (part of which could be caused by the presence of 
the detector) that is not required by a correct quantum 
description of the measurement process. 

It is probably most conventional to describe the qubit 
using a basis where the detector probes the operator cTz 
of the qubit, where c^x,y,z are the Pauli matrices. In other 
words, the detector probes the z component of the pseu- 
dospin associated with the qubit. The qubit Hamiltonian 
then points along an arbitrary direction. For the analysis 
below, we find it easier to describe the qubit in the energy 
eigenbasis. The qubit Hamiltonian is therefore expressed 
as 
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where E is the energy splitting between the qubit's two 
energy levels. We shall express the ground and ex- 
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cited states of the Hamiltonian as |0) and |1), respec- 
tively. The qubit operator that the detector probes 
(i.e. the qubit property that determines the two differ- 
ent switching rates of the detector) can be expressed as 
(Tn = (cos/3)(Tz + (sin/3)(7a;. We shall express the eigen- 
states of this operator as 

\L} = cos^|0) + sin^|l) 

\R) = sin||0)-cos||l). (2) 

We postpone writing down a Hamiltonian of the sys- 
tem because, as we shall see below, the results of Sees. Ill 
and IV are independent of the details of the Hamiltonian 
and the analysis can be performed without knowledge of 
these details, while the behaviour of the detector in the 
case analyzed in Sec. V depends on these details. There- 
fore, by writing down a Hamiltonian at this point we 
run the risk of losing some insight into how the measure- 
ment process proceeds. Some specific Hamiltonians will 
be given in Sec. V. 



A. Simple case 

Let us consider the following simple case for demon- 
stration purposes: We take /3 = 0, such that \L) = |0) 
and \R) — up to an irrelevant phase factor. In this 
case we do not need to worry about any mixing dynamics 
between the states \L) and \R) during the measurement 
process. If we now take the switching rate for one of the 
two states (say \L)) to be exactly zero and the switch- 
ing rate for the other state (here \R)) to be 7, we find 
that one can perform the measurement by connecting 
the qubit to the detector and waiting a duration r that 
is long compared to I/7. If the detector switches to the 
state \S), we know for sure that the qubit was in the 
state \R). If the detector remains in the state |A^), we 
can say with a high degree of c:ertainty that the qubit was 
in the state \L). In other words, instances where 7T :» 1 
and the qubit is in the state \R) but the detector does 
not switch should be very rare and can be ignored for 
all practical purposes. We can therefore reach essentially 
100% measurement fidelity using this procedure. 

The main reason why understanding this case was sim- 
ple is that we did not need to worry about any mixing 
dynamics between the states \L) and \R) during the mea- 
surement process. This simple picture breaks down when 
/3 ^ 0. In that case, one could say that the Hamilto- 
nian tries to push the qubit 's state to precess around 
a certain axis (determined by the qubit 's Hamiltonian), 
while the detector tries to slowly pull the qubit's state 
towards a different axis (determined by the probed oper- 
ator). The interplay between these two mechanisms gives 
rise to nontrivial dynamics and measurement scenarios. 
The analysis of this interplay is the main subject of this 
paper and will be carried out next. 



III. INCOHERENT DETECTOR 

We start with the case where the coherence time of the 
detector (i.e. the dephasing time for coherent superpo- 
sitions of the states |A^) and l^)) is short compared to 
the inverse of the qubit's energy scale E (with the unit- 
conversion factor h) . This limit could also be seen as the 
limit of a classical detector. Note that with the above 
constraint we are not imposing any requirement on the 
relation between the detector's switching rate and the 
energy scale E. 



A. Qubit dynamics during a short period of time 

In order to develop a theoretical understanding of the 
qubit dynamics during the measurement process, we con- 
sider a period of time St that is short enough that (1) 
the detector's switching probability during this period is 
small and (2) the qubit Hamiltonian causes a negligible 
change in the qubit's state, but the same period of time is 
long enough that the switching probability scales linearly 
with time (for time intervals that are short compared 
to the coherence time of the detector, the detector's dy- 
namics resembles coherent oscillations between the states 
I N) and | S) ; for time intervals that are long compared to 
the coherence time of the detector, the detector exhibits 
exponential-decay-like dynamics from the state |A'') to 
the state 15*)). The fact that the detector has a very 
short coherence time is crucial for the simultaneous sat- 
isfiability of the above conditions. Note that, apart from 
its role in deriving some differential equations below, St 
does not have any particular significance. It can therefore 
be chosen freely within the range specified above. 

We identify the short time interval described above by 
its initial and final times, t and t + St. If during this 
interval the qubit is in the state \L), the detector will 
switch with rate 7/,, resulting in a switching probability 
"fLSt. If, on the other hand, the qubit is in the state \R), 
the detector will switch with probability juSt. 

The evolution of the qubit's state can be described 
using the density matrix Pq{t) at time t: 

p,(t)=f^o°w ^oiwy (3) 

We can allow the density matrix to have a trace smaller 
than one (this 'unnormalized' trace would correspond to 
the probability that the dc^tector has not switched from 
the beginning of the measurement until time t) . As men- 
tioned above, we assume that the detector can switch for 
both states \L) and |_R), with rates 71, and (we shall 
generally take 7_r > 7^). 

Assuming that the detector has not switched until time 
t, the state of the entire system composed of the qubit 
and the detector is given by the simple product 

PToUt)=P<,it)(^\N){N\. (4) 
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Because of the qubit-detector interaction, the Hamilto- 
nian takes the density matrix of the entire system into 
the state 

PTotal {t + St) = J/Total {St)pTota\ (t) t^Total i^^) i^) 

at time t + St, where St is taken to be so short that 
all dynamics is coherent (the superscript f indicates the 
transpose conjugate of a matrix; note that several matri- 
ces in this paper will not be unitary, even though some 
of them will be denoted using the letter U). Since the 
time interval under consideration, i.e. t to t + St, is long 
enough that the switching follows exponential-decay-like 
dynamics, we have to modify Eq. ^ such that no co- 



herence between the states \N) and \S) appears on that 
timescale. Making the proper adjustments, and neglect- 
ing for a moment the qubit's free evolution and other 
phase shifts that could be caused by the measurement 
but can be straightforwardly accounted for, we find that 
the evolution of the total density matrix is governed by 
the formula 

PTota.l{t + St) = ?7TotaLs(50PTotal(0^TotaLs(^*) 

+ C/Total,ns(5i)PTotal(0t^Total,„s('^0(6) 

where 



UTot.is{6t) - \L) {L\ ® \S) {N\ + yf^t \R) {R\ ® \S) {N\ 

?7Totai,ns(<5t) = - iLSt \L) {L\ ® \N) {N\ + ^\- jnSt \R) {R\ ® |iV) (iV| . (7) 

One can now imagine that a measurement is performed on the detector at time t + St, projecting it either onto the 
non-switched state \N) or the switched state \S). The state of the qubit evolves accordingly. The qubit's new density 
matrix can be expressed as 

p{t + St)= P,{St)p{t)PHSt) or Pn,{St)p{t)Pl{6t), (8) 

depending on the measured state of the detector (note that we have dropped the index q identifying the qubit's 
density matrix; from now on we only consider the evolution of the qubit's density matrix). The operators Ps{St) and 
Pns{St) can be derived straightforwardly from Eq. ([6]), and we shall give expressions for them shortly. 

To summarize the above argument, the qubit's density matrix evolves in one of two possible ways during the interval 
t to t + St, depending on whether the detector switches or not. The measurement-induced part of the evolution is 
described by the operators Ps{St) and PnsiSt), which correspond to switching and no-switching events during the 
interval of length St, respectively. These projection operators are given by 



PsiSt) = yffiM\L)(L\ + ^^t\R){R\ 

I — J- ( cos^ f sin f cos f \ , / rr ( svc? | - sin f cos | \ 

\ sm I cos f sm | / \ ^ sm f cos | cos"^ f / 



PnsiSt) = - iLSt \L) (L| + v/1 - iRSt \R) {R\ 
2±^i^I^cos/3 ^sin/3 



2 2 ^^'^^^ ,2 

'+1 
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-'-g— sm 1 — cos p 



(9) 



where 7± = (7_r ± 7l) /2. 

The probability that the detector will switch during 
the time interval t to t -\- St va given by 

VxQh^{t,t^St) = i:x\p^{St)p{t)Pl{St)^. (10) 

If the detector switches, the qubit's density matrix im- 
mediately after the switching event will be given by 

p{t + St)^P,{St)p(i)Pl{St), (11) 

up to the usual renormalization constant that is not 
needed for our purposes. The transformation given in 



Eq. pT|) describes a substantial change in the state of 
the qubit. Therefore, the small change induced by the 
qubit's Hamiltonian during this short interval can be ig- 
nored. Note that in this paper we shall not be inter- 
ested in any coherent dynamics that might occur after a 
switching event. As a result, it is not important for our 
purposes that coherence between the states |L) and 
be maintained after a switching event. 

If the detector does not switch between t and t -t- St, 
the qubit's density matrix will be transformed according 
to the operator P^^(St). Since the change induced by 
P-asiSt) is small (proportional to St), one must now take 
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into account the comparably small change induced by 
the qubit Hamiltonian, which is described by the unitary 
transformation 



f/Ham('5t) 



^iE5t/2 







,-iESt/2 



(12) 



The qubit 's density matrix at time t + 5t will therefore 
be given by 

p{t + 5t) « U^U5t)PUmp{t)PL{5t)uLnm- (13) 

Note that tJnamiSt) and Pas{6t) can be treated as com- 
muting operators in the limit 5t — s- 0; as mentioned 
above, we are assuming that 5t is much smaller than 
both l/E and I/7+. 

It is worth pausing for a moment here to comment on 
an issue related to the terminology of strong versus weak 
measurement. There are two possible outcomes of the 
measurement during the interval t to t + St. If the detec- 
tor switches, the qubit undergoes a large change (Eq.fTTj). 
In particular, if 7/, — 0, the measurement would consti- 
tute a strong measurement. A no-switching event, on the 
other hand, causes a very small change in the state of the 
qubit (Eq. [T5|) . Whether such a measurement is weak or 
strong therefore depends on the outcome of the measure- 
ment. One could say that since the strong-measurement 
outcome occurs with a very small probability, the mea- 



surement performed during the interval of length 5t is a 
weak-measurement on average. 



B. Equation of motion for tlie evolution operator 
before switching 

We now take the density matrix of a qubit in an exper- 
imental run where no switching has occurred (yet) and 
express it as 



p{t) = U,,,(t)p{0)Ul(t)- 



(14) 



The operator Uns{t) is therefore the propagator that de- 
scribes how the state of the qubit evolves from the initial 
time t — until time t. Comparing Eqs. and (fH)) . 
we find that 



Uns{t)+ UH.n.{St)Pns{St)-l Uns{t) 



(15) 



Using Eqs. © and (O and taking the limit 6t 0, 
we find that the operator Uns{t) obeys the differential 
equation 



dUnsjt) 

dt 
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This linear differential equation can be solved straightforwardly to give 
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where 

A± 
tan 77 



7+ _^ 1 
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sin/3 



E'^) + 2i£:7_ cos (3 



7_ cos P + iE 



(18) 



Similarly, the qubit 's density matrix after a switching 
event that occurs between t and t + At {At can be un- 
derstood as the accuracy with which the switching time 
can be determined, and for now it is assumed to be much 
smaller than l/E) can be expressed as 



where 



p{t,At) = U,{t,At)p{0)U}{t,At), (19) 



Us{t,^t)^Ps{^t)Uns{t). (20) 



The above expressions will form the basis for the anal- 
ysis of the remainder of this section and that of Sec. IV. 
They describe, in a general setting, how the qubit's state 
evolves during the measurement process. As we shall see 
below, they also describe what information can be ex- 
tracted from a given measurement outcome. 



C. Measurement basis and fidelity for a given 
outcome 

There are a large number of possible outcomes of a 
measurement attempt. One possibility is that the detec- 
tor does not switch throughout a measurement of dura- 
tion r. In this case the corresponding evolution operator 
that gives the qubit's density matrix at time r is simply 
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given by Uns{T) from Eq. ([T7|) . All the other possibilities 
are described by instances where the detector does not 
switch during the interval between the initial time {t = 0) 
and time t and switches between times t and t + At. In 
this case, the relevant evolution operator is Us{t,At). 

The question now is what information (about the ini- 
tial state of the qubit) one can extract from a given out- 
come, e.g. a detector switching event that occurs between 
times t and t + At or a, no-switching instance after time 
r. In other words, what would be the measurement basis 
and fidelity that correspond to a given experimental out- 
come? It should be emphasized here that, in general, the 
measurement basis and fidelity depend not only on the 
measurement procedure, but also on the specific outcome 
obtained in a given experimental run. 

Once a given outcome is observed, one can take the 
relevant evolution matrix Uj, which will be one of the 
matrices Us{t,At) and Uns{t), and think of it as being 
composed of two parts: a measurement operator followed 
by a rotation, 



(21) 



The measurement operator takes the form 



7- rmcas 
j 



(1) 



where 



(1) 



p. 
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(2) 



0} 



(2) 



(22) 

are two orthogonal states that 



represent the measurement basis of this particular out- 
come, and Pj^^ and Pj"^^ are, respectively, the probabil- 
ities that the specific outcome under consideration will 

be observed if the qubit was initially in state i/i^^^ 

;(2) 

For a given initial state p{0), the probability that the 
specific outcome under consideration occurs in an exper- 
imental run is given by 



Prob, = 



Tr 
Tr 

Tr 



^ P(0) 



^mcas I 



(23) 



the maximally mixed state 



Pinax. mixed ( Q 2. 



If we choose the labels in Eq. ^ such that pj^' > P^ 
one would maximize the probability of making a correct 
inference about the qubit 's state by associating the spe- 



(24) 
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(i; 



Using the max- 



cific outcome j with the state 

imally mixed state, a straightforward calculation shows 
that the fidelity in this case is given by 



P, 
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(25) 



When there are more than one possible outcome for the 
measuring device (which is always the case for any useful 
measuring device), the overall fidelity can be evaluated 
by taking the statistical average of the fidelities for the 
different possible outcomes. This averaging procedure 
will be performed in Sec. IV. We should emphasize, how- 
ever, that the fidelities of the different specific outcomes 
also carry some significance of their own. This point will 
become clearer when we discuss a number of illustrative 
examples below. 

An alternative approach to deriving the fidelity goes 
as follows: let us assume that we initially have no infor- 
mation about the state of the qubit. The state is there- 
fore described by the maximally mixed state in Eq. (|24p . 
After the detector gives the outcome denoted by the in- 
dex J, we have partial or complete knowledge about the 
qubit 's state immediately after the measurement, or in 
other words the collapsed qubit state. The amount of in- 
formation that we have gained can be quantified by the 
purity of the qubit 's state after the measurement; 



Purity = W2Tr 



max. mixei 



- 1. 



(26) 



A straightforward calculation shows that this expression 
agrees with the one derived above for the measurement 
fidelity. 



As a result, by diagonalizing the measurement operator 
jjmoas ^Qj, equivalently, by diagonalizing ujUj), we can 
obtain the measurement basis and fidelity of the per- 
formed measurement. 



The measurement basis is defined by the states 







(1) 



and 



The measurement fidelity is given by the 



difference between the probability of making a correct 
inference about the state of the qubit and the probability 
of making a wrong inference about the state of the qubit. 
In order to describe all possible input states, one must 
calculate the fidelity by taking a statistical average using 



IV. ANALYZING DIFFERENT IMPORTANT 
CASES 

The general expressions that we have derived in Sec. Ill 
do not provide a clear, intuitive picture of how the qubit 's 
state evolves during the measurement and how the mea- 
surement data should be interpreted. We therefore an- 
alyze some special, representative cases below. Section 
IV. A analyzes the simple case discussed in Sec. II. A, 
while the choice of the cases analyzed in Sees. IV. B and 
IV. C is prompted by the fact that Eq. ^TE\\ can be sim- 
plified in these opposite limits. 
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A. Case 1: P — (No mixing between the states \L) 
and \R) during the measurement) 

The case /? = is the case where the Hamihonian 
and measurement axes are paraUel (in other words, the 
states \L) and \R) coincide with the states |0) and |1)). 
We can therefore expect to reproduce the simple, known 
resuhs outhned in Sec. II. A. When /3 = (which also 
gives 77 = 0), we find that 
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Wc therefore find that 

Unsit) = 
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]Ur''%t,At), (29) 



where 



L/— (i) 



e 2 




VlT^te-^' 
VlR^e 



-2^t (-30) 



The eigenvectors of the matrices U^^'^'^^{t, At) and 
U^'^'^^^t) are the states \L) and \R) (independently of the 
values of t and At). The measurement basis is therefore 
\R)}, as expected. The eigenvalues of U^'=''''{t, At) 
are given by 



pL{t,t + At) = VT^Atcxp 



'P,Rit,t + At) = y/jRAtCXp 
and those of U^°'^^{t) are given by 



2 

lRt\ 

' 2 i 



(31) 



PL{t) = exp 



pR{t) = exp 



2 J 
2 / 



(32) 



Fidelity of a specific outcome 



The eigenvalues in Eq. (pij) give a fidelity (for switching 
time t) of 



Fsit,At) 



P,^{t,t + At)-P,Rit,t + At) 



pL{t,t + At) + pR{t,t + At) 



-iLt , 



-fB.t 



(33) 



If the detector switches soon after the measurement 
starts, i.e. much earlier than tq = (log7fl — log7L)/(7ii; — 
7l) [with log being the natural logarithm], the likely state 
of the qubit is \R) with fidelity 



F,(0,At) 



IR - IL 



(34) 



whereas if the detector switches at a time much later 
than Tq, the likely state of the qubit is \L) with a fi- 
delity that approaches one. Interestingly, if the detector 
switches at time tq, the fidelity vanishes. In this case 
the measurement does not give any information about 
the state of the qubit, and the state of the qubit immedi- 
ately after the switching is equal to the initial state up to 
a simple rotation [lol |. In other words, the purity of the 
qubit 's state does not change when this outcome is ob- 
served (note that the qubit's state is, in general, a mixed 
state). The function Fs{t, At) is plotted in Fig. 2 for the 
case 7_r/7l = 10. 




FIG. 2: The measurement fidelity Fs(f, At) as a function of 
the switching time t for 'ypi/'yL = 10 and /3 = 0. The accuracy 
to which the switching time is known, At, is assumed to be 
much smaller than 1/E. Therefore, the fidelity is independent 
of At and is a function of t only. Recall that to = (log7ij — 
log7i)/(7fl - 7i)- 



2. Overall fidelity 

The overall fidelity is now evaluated by taking the sta- 
tistical average over all possible outcomes 



F = 



^ ^ Tl t^j Pmax. mixed 

j 

F.it^dt)&2l±^^^^±&2l±^ 



iRe - 7Le 



8 



(7fle-^«* + 7Le-^^*) dt 



(7Le-^^* - 7fle-^«*) 



To 



e " — e 

7-R \ "'-"""'^ I IR^ -yR-iL 



IL 



(35) 



This expression for the fidehty is plotted as a function 
of IrHl in Fig- 3. One could say that for large values 
of IrIil one has a strong measurement, whereas small 
values of ^rJ^l correspond to weak measurement. As 
can be seen from Fig. 2, however, even for small val- 
ues of 7fl/7L there can be instances where a high-fidelity 
outcome is obtained (when t/ro 3> 1). 




40 60 

IbJil 



100 



FIG. 3: Overall fidelity as a function of the ratio of switching 
rates ')r/^l when the qubit Hamiltonian commutes with the 
probed qubit operator, i.e. in the special case /3 = 0. 

It is interesting to compare the above result for the 
overall fidelity with the one that would be obtained in 
the following situation: Let us assume that we connect 
the qubit to the detector for some duration r and only 
at the end check whether the detector has switched or 
not. A switching instance would be associated with the 
state \R) and a no-switching instance would be associated 
with the state \L). The relevant switching probabilities 
are given by: 



Psi^^r) - (1 
^s^(0,t) ^ (1 



-lL'r\ 



-lRr\ 



(36) 



The fidelities for switching and no-switching instances 
are therefore given by 



= 1LT 



2 - e-Ti^ - e-iR'' 



(37) 



and the overall fidelity is given by 



F = 



-IRT 



2 - e^'^^'' - e^T-"^ 



-IRJ 



= tl''^ _ a IR'^ 



2 



(38) 



This expression for the overall fidelity is maximized by 
choosing the pulse duration r = tq (which was defined in 
Sec. IV.A.l). This means that the optimal pulse dura- 
tion is determined by the switching time tq at which the 
fidelity Fs(0, tq) of Eq. p4)) vanishes, i.e. the switching 
time that separates the two time regions associated with 
the states \L) and \R). No enhancement in overall fidelity 
is obtained by keeping track of the exact switching time. 

One thing that can be gained by keeping track of the 
exact switching time is our knowledge about the state 
of the qubit after the switching event in a specific ex- 
perimental run. For example, starting with a maximally 
mixed state and choosing instances associated with high 
fidelity (e.g. those with no switching until times much 
later than tq), one can predict with a high degree of cer- 
tainty the qubit's state after the measurement. 



B. Case 2: S <C 7- (Instantaneous measurement) 

This case is quite simple as it corresponds to the 
limit of instantaneous measurement [note that the above 
condition guarantees that E <C ^r] recall that 7_ = 
{iR — 7l)/2]. Equations now reduce to 



7+ , 1 

± -7_ 

2 2 ' 



(39) 



such that 



Usit,At) 



rrmcas 
^mcas 



it) 

{t, At) 
\L) {L\ 



\R) {R\ 



yJlhAt 



-)-/^e-**|i?)(i?|. (40) 



Note that the energy E does not appear in the above 
expressions; the detector switches so fast, at least for one 
of the qubit states, that the measurement is completed 
before the qubit Hamiltonian contributes any change to 
the state of the qubit. 

The fact that E can be ignored in this case is simi- 
lar to the situation analyzed in Sec. IV. A (case 1). One 
can also see that the last two lines in Eq. (|40)) coincide 
with Eq. ([50]) . One can therefore follow the analysis of 
Sec. IV. A [starting from Eq. (|30p ] in order to calculate the 
measurement fidelity for the instantaneous-measurement 
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case under the different possible scenarios. Not surpris- 
ingly, the measurement basis is {\L) , which is natu- 
ral given the fact that the measurement is instantaneous, 
and the measurement basis should coincide with the nat- 
ural measurement basis of the detector. 

As in Sec. IV. A, the overall fidelity is given by Eq. ([35]) 
[which is plotted in Fig. 3]. In the limit ^r/^l oo, the 
fidelity is one, and one obtains the simple case of a strong 
(projective) measurement in the basis {\L) , Note 
that the limit Jr/^l gives the same result, with the 
roles of and jr reversed. 

In the hmit 7_ <^ 7+ (i.e. when ^r/^l is close to one), 
we still have an instantaneous measurement, but with a 
low fidelity. The measurement in this limit is therefore a 
weak measurement. The switching time of the detector 
allows the experimenter to make a guess about the state 
of the qubit, but with a small degree of confidence. Re- 



peated measurements of this type, possibly synchronized 
with the qubit 's free precession, on a single qubit would 
increase the acquired amount of information (or degree of 
certainty) about the state of the qubit, thus approaching 
the case of a strong measurement. 



C. Case 3: £ 3> 7+ (Slow measurement) 

The opposite limit of case 2 is the one with ^ 7_ . 
However, when 7_ <C ^ 7+, we recover the case of an 
instantaneous, weak measurement (with minor changes). 
Therefore, in this subsection we focus on the case where 
E is larger than both 7_ and 7+. 

When £' 3> 7- , we find that (see Appendix A) 



uLmnsit) 



g-(7+-7- cos/3)t 







-(7++7- cos/3)t 



(41) 



where 



e(t) _iEf 
cos -j^e 2 

■ Sit) iE-t 

sm -^e 2 



sm -^e 



cos -j^e 2 



E = E-jl/2E 



e*-' 



e't< 



cos 
sin 



m 

eft) 



m 



tane'(t) 

pW(i,Ai)+Pi2Hi,Ai) 



{iL - 7fl)sin/3 



7L (^eT-'^^/^cos^ I - e-'>'-*™^'3sin^ |^ + 7^ (^eT-*™'^/' sin^ | - e-^-*™'''^ cos2 |^ 
7L (e^-'-^'f'cos' ^ -f e-^-'^^^'sin^ ^ +^r (e''-"'°'l'sm^ ^ + e"^-*™'''' cos^ ^ 



ilL - 1r) sin/3 
sin6'(t) 



(42) 



Recall that, in spite of their appearance, the operators 
Uns{t) and Us{t,At) are not unitary matrices. 

For short times {t — 0), the above expressions reduce 

to 

e ^ p 

pW+Pi^) = Ai(7L+7fl) 

p(l) _ p(2) ^ ^^(^^ _ 



P,(0,At) = , '\ (43) 

such that the measurement basis coincides with the de- 
tector's natural measurement basis, as it should for short 
times (where the qubit Hamiltonian has not affected the 
qubit's state). For long times {j-t 00) [and (3 ^ 7r/2], 
on the other hand, the above expressions reduce to 
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0^0 

Pi^) + = e-^+*A< |7Le^-*™^'' cos^ | + 7^6^-*™^'' sin^ | 

F{t'^oo,At) = 1, (44) 



such that the measurement is performed in the qubit's 
energy eigenbasis with 100% fidehty. Thus the measure- 
ment basis starts from the detector's natural measure- 
ment basis (i.e. {\L) , \R)}) at t = and spirals towards 
the z axis (i.e. {|0) , |1)}) with time. Note that since the 
denominator in the first line of Eq. (|42p goes from neg- 
ative values (e.g. at t = 0) to positive values (e.g. when 
t 00) the measurement basis goes through the equator 
during its spiraling motion. 

From the above results, we can see that the measure- 
ment is made in a rotated basis that is determined by 
the exact switching time, which is uncontrollable. Each 
experimental run gives a measurement outcome that is 
different from other identically prepared runs. Note, in 
particular, that when 7l = the fidelity is 100% for all 
switching times t (which can be verified with straight- 
forward algebra). This result is quite interesting; even 
though the fidelity is 100% for every single run, the mea- 
surement basis is unpredictable and is only determined 
when the switching event occurs. 

Let us now consider the situation where we know only 
whether the detector switched between t = and t = t 
(with Et ^ 1) or not. The measurement matrices that 
correspond to these two possible outcomes are given by 



U^,{0,t)Us{0,t) = / ijl{t,dt)U,{t,dt) 



= l-Ul{T)UUr). 



(45) 



and i/2s(r)J7ns(''') as given in Eq. (|¥T|) . Since the measure- 
ment matrices are diagonal in the qubit's energy eigen- 
basis, the measurement is performed in that basis. The 
overall fidelity can be calculated similarly to what was 
done in Sec. IV. A: 

p — g-(7+-7- cos/3)t _ g-(7+-|-7_ cos/3)t ^-^g-j 

This expression is maximized by choosing 

o -y_|_ — COS p 

27„ cos (3 

which gives the maximum fidelity 



(47) 



7+ -7- 
7+ — 7_ cos (i \ '^'<- 

7+ + 7- cos /3 

7+ + T_ 

( 1+ ~1~ cos /? ^ 
V 7+ + 7- cos (3 



• (48) 



Since the above expression looks somewhat cumbersome, 
we shall not pursue it in its general form any further. In 
the special case 7^ = 0, the overall fidelity for a readout 
pulse of duration r is given by 



p — g-7-R-rsin ^ _ g-7fircos ^ 



which is maximized by choosing 

— log tan^ I 
7i?^ cos (3 

The maximum fidelity is therefore given by 



(49) 



(50) 



2sin^(l3/2) 2coB^(/3/2) 



tan 



/3 



sec (3— 1 



tan 



sec /3 -1-1 



(51) 



This expression is plotted in Fig. 4. The measurement 
fidelity decreases as (3 approaches 7r/2 and vanishes at 
that point. 




FIG. 4: Maximum fidelity as a function of the angle between 
the qubit Hamiltonian and measurement axes, assuming that 
the exact switching time is not known (up to the measurement 
duration, which is taken to be much larger than l/E). 

For further demonstration of the results discussed in 
this section, see Appendix B, where the special case (3 — 
7r/2 is analyzed in some detail. 
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D. All-in-one measurement setting 

Looking at the results of Sec. IV. C (case 3), it is in- 
teresting that every time the experiment is repeated, the 
measurement outcome gives information about the qubit 
state in a different basis (that cannot be predicted be- 
forehand). This result suggests that it could be possi- 
ble to perform 'single-setting quantum state tomography' 
[ni; instead of the externally apphed rotations used to 
perform measurements in different bases in conventional 
quantum state tomography, the (uncontrollable) switch- 
ing times determine the measurement basis in each ex- 
perimental run. Indeed, using Eq. (|64p we find that the 
switching probability as a function of time is given by 

Prob.(<,A<) = Tr{f7,^(i,Ai)f7s(t,At)p(0)} 

= e-^+*At pooe^~*™^''(7+-7-Cos/3) 
+pne-^-*^°^'^(7++7- cos/3) 



7— 

(/Oio + pai) —sm(3 cos Et 

( PlO - POl \ 7- . r, ■ ^, 

— smpsnii/i 



(52) 



By fitting the experimentally measurable switching-time 
probability distribution to the above function, one can 
extract the three quantities (pn — poo), (pio + Poi) and 
(pio ~ Poi) /i^ which completely characterize the quan- 
tum state of the qubit. In addition to the above three 
parameters, one can also extract the values of 7l, 7_r, 
(3 and thus providing an all-in-one measurement that 
characterizes both the systern parameters as well as the 
quantum state of the qubit [I4I. Note that in the spe- 
cial cases (3 = 0, /3 — 7r/2 and 7+ ^ full system 
characterization and state tomography are not possible, 
in agreement with intuitive expectations. In particular, 
the case [3 = Q corresponds to a measurement in the 
{|0) , |1)} basis (regardless of the switching time), which 
clearly cannot give any information about the quantities 
(pio + Poi) and (pio - poi) /«■ 



V. COHERENT DETECTOR 

In the preceding two sections, we have focused on the 
regime where the coherence time between the two detec- 
tor states (i.e. the switched and non-switched states \S) 
and |A^)) is the shortest timescale in the problem. Now 
we turn to the case where this time is longer than 1/-E, 
which now is the shortest timescale in the problem. 

We should emphasize here the distinction between this 
case of a coherent detector and the case of a 'slow' but in- 
coherent detector (with a small switching rate; 7+ ^ E) 
analyzed in the previous section [l^. That situation 
was analyzed using the same theoretical framework as 
any other case with an incoherent detector. We should 



also emphasize that, in spite of the detector's long coher- 
ence time (compared to i/E), we still assume irreversible 
switching. These assumptions are compatible; for irre- 
versibility one only needs to assume that the detector's 
coherence time is much shorter than the inverse of the 
switching rate. The assumptions explained above can be 
summarized as follows: 

— ^ coh. time of detector ^ — — — . (53) 

E switching rate 

Since we are now treating the detector as a quantum- 
coherent object (at least on the timescale of qubit preces- 
sion) , we analyze the dynamics using a Hamiltonian that 
describes the combined qubit-detector system. Since the 
Hamiltonian contains the physical interaction between 
the qubit and the detector, in this section we find it eas- 
ier to define the qubit's z axis as defined by the probed 
operator. Thus, the Hamiltonian can be expressed as 



dct •) 



(54) 



where 



-— (sin/3(7j; + cos/Jo-^) 

{gL{\S) {N\ + h.c.)+eL\S) {S\}\L) {L\ + 
{gR {\S) {N\ + h.c.) + en \S) {S\} \R) {R\ ,(55) 



gL and g^j are the coupling strengths that induce detector 
switching for the qubit states \L) and respectively, 
and and are relative energy biases between the 
states |A^) and \S) for the two qubit states. As before, 
the detector is assumed to be initially in the state \N). 

Since the qubit's energy is the largest energy scale 
in the problem, the dynamics described by the above 
Hamiltonian and the detector's decoherence cannot in- 
duce transitions between the qubit's ground and excited 
states (note also that we are ignoring any intrinsic de- 
coherence mechanisms of the qubit). We can therefore 
analyze the dynamics by treating two separate subspaces 
of the Hilbert space, one with the qubit in the ground 
state |0) and one with the qubit in the excited state |1). 
This approach results in two effective Hamiltonians for 
the detector: 

^^dct,|0) = qubit(0| |0)q^bit 

= -go{\S){N\ + h.c.) + lo\S){S\ 

■f^dct.ll) = qubit(l|-H"|l)qubit 

= -gi{\S){N\ + h.c.) + -ei\S){S\, (56) 



where 











/3 




= 5lcos^ 
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f gR siv? 


2 






P 




/3 


91 


= 5Lsin^ 
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V gn cos^ 


2 






P 




/? 


eo 


= e^cos^ 


2 


1- eji siv? 


2 










f3 


ei 


= €l sin^ 




- eR cos^ 


2 



(57) 



In order to avoid running into complicated expressions 
below, we make the experimentally relevant assumption 
that the switching for the state \R) is much faster than 
that for the state \L). This situation occurs when one or 
both of the following conditions are satisfied: (1) \gR\ ^ 
I^lI or (2) jeil is much larger than all other parameters 
relevant to the switching dynamics. We analyze these 
two cases separately. 

First we consider the case where \gii\ ^ \gL\, such 
that gL can be ignored in the following analysis. For 
simplicity and clarity, we take and ej^ to be smaller 
than the detector's coherence time, such that they also 
can be ignored for purposes of this calculation. Equations 
([56]) therefore reduce to 

^det.io) = gB^sin^^{\S){N\ + h.c.) 

^dct,|i> = 9R cos2 ^ (1^) (N\ + h.c.) . (58) 

Since we are assuming that the switching is an incoher- 
ent process (or, in other words, the detector's coherence 
time is much shorter than |(7i?|), the switching rates are 
proportional to the squares of the transition coefficients 
in the above Hamiltonians. The switching rates for the 
qubit states |0) and |1) can therefore be expressed as 

-n -n ■ 4 /3 

Fo = Fi^sm - 

Fi = Fficos^^. (59) 

We now consider the second case, where le^l is very 
large compared to all other parameters relevant to the 
switching dynamics. In this case the switching rates are 
governed by the scaling formula 

With the simplifying assumption that = gn^ we find 
switching rates of the form 

Fo = Fflsec^l 

Fi = Fflcsc^^. (61) 

Note that the divergence of Fi for /3 is an artefact 
of our approximations. In a more accurate description, 
Fi would saturate at a finite value as [3 approaches zero. 

The difference between the expressions in Eqs. (j59p 
and (|6T|) shows that, in the case of a coherent detector, 
knowledge of the microscopic details of the detector is 
necessary in order to determine the dependence of the 
switching rates on the angle (5. On the other hand, this 
sensitivity to details could be useful for experimentally 
testing theoretical models of the detector. 

One common, and important, result for the different 
cases with a coherent detector, however, can be seen 
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Detector type 


Measurement 
basis 


Meas. fidelity 
limited by 


Fast switching 


{\L)AR)} 




Slow switching, 
incoherent, switching time 
accurately measurable 


determined 
upon switching 


Ir/il 


Slow switching, 
incoherent, exact switching 
time inaccessible 


{|0),|1)} 


7h/7l and /3 


Slow switching, 
coherent 


{io),|i)} 


7^/71- and (3 
(in two-state 
model; Sec. V) 



Table I: Measurement basis and fidelity for diflferent types of 

detectors characterized by their coupling strength to the 
qubit, their coherence times and whether the switching time 
can be known accurately or not. In the case of a weakly 
coupled coherent detector, we give the results for the 
two-state-detector model only. 

by noting that the qubit 's state in the energy eigenba- 
sis {|0) , |1)} does not change during the measurement. 
Therefore, the detector cannot be performing a measure- 
ment on the qubit in any basis other than this one. Since, 
in addition, the detector's switching rate depends on the 
qubit's state in the energy eigenbasis, one can conclude 
that the detector does in fact perform a measurement on 
the qubit in the energy eigenbasis. Using the transition 
rates given in Eq. ([59]) or ((6T|) , one can follow the analysis 
of Sec. IV.A and determine the fidelity of the measure- 
ment setup. The two expressions found in this section 
describe a fidelity that slowly decreases from unity at 
/3 = to zero at /3 = 7r/2. This result seems to be un- 
avoidable for the two-state-detector model considered so 
far. In spite of the different functional dependence, the 
fidelity decrease is somewhat similar to the one obtained 
in Sec. IV. C for an incoherent detector with small switch- 
ing rates [see Eq. (fSTjl]. In Sec. VII, we shall discuss a 
system where the fidelity of a coherent detector remains 
very close to unity even as (3 approaches 7r/2. 

VI. COMPARISON BETWEEN DIFFERENT 
TYPES OF DETECTORS 

To conclude the above analysis of the detector's oper- 
ation, we present in Table I the measurement basis and 
the factors limiting the measurement fidelity in different 
operation regimes. 



It is worth mentioning here a distinction between two 
different types of weak measurement, i.e. a measurement 
with low fidelity. There are cases where the interaction 
between the qubit and the detector is such that it is fun- 
damentally impossible for any observer (with full access 
to all degrees of freedom of the detector) to determine 



13 



with certainty the state of the qubit. The situations ana- 
lyzed in Sees. IV. A, IV. B and V correspond to this type of 
weak measurement. In this case, there is no lost informa- 
tion, and an initially pure state of the qubit remains pure 
after the measurement (note that not knowing the exact 
switching time to an accuracy that is much smaller 
than 1/ E does not constitute lost, and relevant, informa- 
tion in this context; Going to higher accuracies within 
the regime ElSi <C 1 provides more information about 
the state of the detector, but not that of the qubit). The 
situation where the observer cannot establish the exact 
switching time of an incoherent detector to an accuracy 
smaller than 1/E (analyzed in Sec. IV. C) involves infor- 
mation that is collected by the detector and potentially 
measurable, but inaccessible to the observer in the exper- 
imental situation under consideration. In this case, the 
qubit 's state after the measurement is generally a mixed 
state (obtained by averaging over all possible outcomes 
within the experimenter's accuracy range), even if the 
initial state is pure. One could say that the accessible 
information serves as measurement data, whereas infor- 
mation that is collected by the detector but is inaccessi- 
ble to the experimenter plays the role of a decoherence 
mechanism. 

We now pause for a moment to comment on the is- 
sue of qubit decoherence. In some experimental situ- 
ations, the detector enhances the qubit's decoherence 
rate [ij]. This additional decoherence channel typically 
causes loss of coherence in the {\L) , basis. If the 
decoherence rate is larger than the qubit's energy scale 
i?, the qubit's state will be frozen in the {\L) , |_R)} ba- 
sis during the measurement. In that case, the detector's 
operation resembles that in the instantaneous (and pos- 
sibly strong-measurement) regime, even if the switching 
rates are small compared to E. Thus, such an additional 
decoherence could be detected by analyzing the opera- 
tion of the detector and comparing it with the expected 
behaviour in the absence of decoherence [lB|. It is also 
worth mentioning here that including decoherence into 
our theoretical analysis would be quite nontrivial. The 
evolution of the qubit's state during the short interval 5t 
would have to be described by the sum of terms, and not 
just a single propagator as was the case in our analysis. 
It would therefore not be possible to derive a simple an- 
alytic expression for the matrix U-as{t) in the presence of 
decoherence as we were able to do in Eq. (|17l) . 

To conclude the discussion of the two-state model, we 
summarize some of the main principles that have been 
demonstrated in the above analysis. For a strong mea- 
surement, the detector always measured the qubit in the 
, |i?)} basis. In the case of a slow, incoherent detec- 
tor whose switching time can be determined accurately, 
one could say that, to a lowest approximation, the qubit's 
state precesses freely following the qubit Hamiltonian and 
is then measured in the {|L) , basis when the de- 
tector switches. When the switching time can only be 
known with small accuracy (i.e. much longer than the 
qubit's precession period), the detector effectively probes 



the qubit in the energy eigenbasis, with the correspond- 
ing switching rates being obtained as weighted averages 
of the bare switching rates 7l and ^r. A slow, coher- 
ent detector measures the qubit in the energy eigenbasis, 
with the detector effectively following weighted-average 
Hamiltonians for the two energy eigenstates. 



VII. APPLICATION TO SUPERCONDUCTING 
FLUX QUBITS 

In the preceding sections we have analyzed a theo- 
retical model where the detector is a two-state system. 
We now turn to a specific physical example of switching- 
based measurement, namely the readout of a supercon- 
ducting flux qubit [l6|. Although the two-state-detector 
model does not provide an accurate description of the 
complex circuits that we discuss in this section, the in- 
sight developed in the preceding sections (in particular 
the general principles summarized in Sec. VI) will allow 
us to infer a number of results concerning the realistic 
situations that we consider here. 

There are two qubit-readout methods used in exper- 
iment that rely on a switching process in a supercon- 
ducting circuit: switching from a zero- voltage to a finite- 
voltage state across a dc superconducting quantum inter- 
ference device (SQUID) [13| and switching between two 
possible dynamical states of a nonlinear oscillator [l^ 
(For a theoretical stud y o n measurement in supercon- 
ducting qubits, see e.g. [l9|; and for studies on quantum 
tunneling between dynamical states in superconducting 
systems, see e.g. [lo. Hi], H^). Both of the readout 
methods mentioned above can be modelled theoretically 
as a fictitious particle moving in a trapping potential. 
Therefore, their behaviour in terms of switching dynam- 
ics should exhibit the same features, which we discuss 
below. 

We briefly explain the operation of the dc-SQUID de- 
tector that relies on switching from the zero- voltage state 
to the finite-voltage state of a current biased dc-SQUID 
(see also [5ll23|); as mentioned above, the analysis of the 
dynamical bifurcation would be similar. The principle of 
the measurement process is as follows: The dc-SQUID 
is known to be a highly sensitive device for measuring 
magnetic fields. In the presence of an externally applied 
magnetic flux and a bias current that is close to a cer- 
tain (critical) value, the switching rate of the dc-SQUID 
from the zero-voltage to the finite-voltage state varies 
greatly with even small changes in magnetic field. When 
a flux qubit is placed next to a dc-SQUID biased close to 
the switching point, the clockwise- and counterclockwise- 
current states of the qubit can result in two switching 
rates that are different by orders of magnitude. By choos- 
ing suitable values for the applied bias current and the 
length of the current pulse, the SQUID switches with 
high probability for one qubit state and with low prob- 
ability for the other qubit state. The contrast between 
the two switching probabilities is perhaps most clearly 
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visualized in the so-called S curves. These curves are 
produced by plotting the SQUID 's switching probability 
as a function of bias current for the qubit's clockwise- and 
counterclockwise-current states (for a fixed measurement 
duration) . An example of such S curves is shown in Fig. 5 
(dashed hues). 

We now ask the question of what the S curves would 
look like if the qubit were biased close to the so-called de- 
generacy point, where the energy eigenstates are super- 
positions of the two persistent-current eigenstates. Fig- 
ure 5 shows three possibilities that correspond to three 
different types of detectors: a strongly coupled detector 
(see Sec. IV. B), a weakly coupled fast-decohering detec- 
tor (see Sec. IV. C) and a weakly coupled slow-decohering 
detector (see Sec.V) [the explanation of the reasoning be- 
hind each one the obtained shapes of the S-curves is given 
in Sec. VI]. In Fig. 6 we plot the measurement fidelity (at 
the optimal bias current) for the three different types of 
detectors. Clearly, these three types of detectors operate 
drastically differently. The fidelity in Fig. 6(a) is given 
by cos/3, and the one in Fig. 6(b) is given by Eq. ([5T|l . 
i.e. it is identical to the one shown in Fig. 4. As men- 
tioned above (see also Q), the measurement fidelity of a 
coherent detector (Fig. 6c) can remain close to 100% even 
when the qubit Hamiltonian and the probed operator do 
not commute. 

The operation of the detector when the qubit is bi- 
ased close to the degeneracy point can be used to identify 
whether it performs a strong- or weak-coupling measure- 
ment and how the detector's coherence time compares to 
the qubit's precession period. By tuning the qubit's en- 
ergy splitting and probing it using the detector, it could 
be possible to measure the coherence time of the detec- 
tor. This approach could also be used to probe coherent 
behaviour in the detector. 

Finally, we comment on the description of a 
continuous- variable switching detector, which should be 
a good description for the dc-SQUID-based measurement 
of the flux qubit (see also Q). Instead of our two- 
state-detector model, one can in general use a model 
of a fictitious particle initialized in a local minimum 
of a one-dimensional potential. The clockwise- and 
counterclockwise-current states of the qubit correspond 
to slightly different potentials for this particle. The two 
potentials can result in very different escape rates (out 
of the local minimum). We expect that in the case of 
an incoherent detector the above results for a two-state 
detector also hold for the continuous-variable detector. 
In the case of a coherent detector, one can follow the 
general procedure that we used for the two-state detec- 
tor, modified for the continuous-variable detector. One 
takes the two trapping potentials for the clockwise- and 
counterclockwise-current states. From these potentials 
one derives two effective potentials (that correspond to 
the qubit's two energy eigenstates) for the fictitious par- 
ticle, each potential obtained as a weighted average of 
the two bare potentials. The proper proportions of the 
original potentials that are used to evaluate the energy- 



eigenstate effective potentials are determined by the mix- 
ing probabilities of the persistent-current states in the 
energy eigenstates. These weighted-average potentials 
now determine the relevant escape rates for the parti- 
cle, and these escape rates characterize the operation of 
the detector. As in Sec. V (two-state coherent detec- 
tor), the detector performs a measurement on the qubit 
in the energy eigenbasis. However, as can be seen in 
Fig. 6(c), the continuous- variable detector has the ad- 
vantage that it could lead to a high fidelity measurement 
even when the energy eigenstates are superpositions of 
persistent-current states. In fact, such high-fidelity read- 
out in the energy eigenbasis was observed in the experi- 
ment of Ref. Q . 

A note is in order here on the relevant coherence 
time that separates the coherent-detector and incoherent- 
detector regimes for a continuous- variable detector. Even 
if the coherence time between the switched (or escaped) 
state and the non-switched (or trapped) state is short, 
the coherence time that corresponds to motion around 
the local minimum of the trapping potential could still 
be long. If (1) the temperature is low enough such that 
the detector's fictitious particle is almost in the ground 
state of its effective potential and (2) the coherence time 
for superpositions of the ground states of the two possi- 
ble effective potentials is longer than 1/i?, the detector 
would operate in the coherent regime [Fig. 6(c)]. A full 
understanding of the boundary at which decoherence in 
the detector causes irreversible collapse of the state of 
the system and how that affects the state of the qubit 
is lacking. Our results demonstrate that the coherence 
properties of the detector are crucial for determining its 
regime of operation. However, they do not answer ques- 
tions about the coherence within the trapped state, or 
states, of the detector. 



VIII. CONCLUSION 

We have analyzed the operation of a switching-based 
detector that is used to measure the state of a qubit in 
the case when the probed observable does not commute 
with the qubit Hamiltonian. We have found several op- 
eration regimes depending on the relation between the 
detector's switching rate and the qubit's energy (strong 
versus weak-coupling measurement) and the relation be- 
tween the detector's coherence time and the qubit's en- 
ergy (coherent versus incoherent detector). The acces- 
sibility of the exact switching time also plays a role in 
determining the performed measurement. Table I sum- 
marizes the differences between the different regimes. 

The weak-coupling regimes result in a number of inter- 
esting phenomena. For example, an incoherent detector 
can be used for an 'all-in-one' measurement setting that 
can be used to simultaneously calibrate the system and 
perform quantum-state tomography. Apart from isolated 
special coherent detector always measures the 

qubit in its energy eigenbasis, regardless of the quantity 
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that it naturally probes. Furthermore, the measurement 
fidelity of a coherent detector can remain close to 100%, 
even when the qubit's Hamiltonian and the qubit opera- 
tor being probed by the detector do not commute. 

Although we have analyzed a simplified two-state- 
detector model that was inspired by the measurement of 
superconducting flux qubits using dc SQUIDs, the main 
results should apply to other types of qubits and mea- 
surement devices, e.g. charge qubits measured using a 
quantum point contact [13] . As far as the roles of the dif- 
ferent coherence and coupling parameters are concerned, 
our results provide insight into the operation of a gen- 
eral type of detector and can be useful for future studies 
aimed at designing high-fidelity measurement devices. 
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In order to determine the measurement basis for instances where a switching event occurs between t and t + At, 
the matrix {t, At)Us{t, At) must be diagonahzed, i.e. converted into the form 
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Clearly (p = Et, such that the measurement basis rotates about the z axis with frequency E, similarly to the free 
evolution of qubit states but with a slightly reduced rate (note also that the sense of rotation is opposite to that of 
state precession). With straightforward algebra we find that 
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Appendix B: Equation P7|) therefore gives 

Incoherent detector with 7+ and /3 = 7r/2 

For further demonstration purposes, we analyze in this 
Appendix the case of an incoherent detector with E 3> 
7+ and /3 = 7r/2. In this case Eq. pB]) reduces to 
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where i? = y — 7^ . The above two matrices can now 
be re-expressed as 
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The interpretation of these results is obvious. In the 
case of no switching, the qubit state evolves freely accord- 
ing to the Hamiltonian (i.e. rotates around the z axis by 
an angle approximately given by Et), but we essentially 
do not gain any information about the state of the qubit. 
On the other hand, an experimental run with the out- 
come of no switching until time t followed by a switching 
event between times t and t + At is equivalent to a mea- 



surement in the basis 



(t)) {t))} followed by 



a rotation around the z axis by an angle Et. The fidelity 
of the measurement is given by 



F,{t,At) 



(71) 



It is interesting that the fidelity is constant in time. The 
only change that occurs in time is that the measurement 
is performed along a rotated axis that precesses about the 
axis of the qubit Hamiltonian at rate E. As explained 
above, the measurement axis is determined only at the 
time that the detector switches. Of course an equivalent 



description of the measurement dynamics would be to say 
that the qubit state precesses freely around the z axis by 
an angle Et and is then measured in the , basis 
when the detector switches. Depending on whether one 
is interested in reaching an intuitive description of the 
qubit-state evolution or in identifying what information 
can be extracted regarding the initial state of the qubit, 
one or the other of the above descriptions would be more 
relevant. 



We can now use the above results to demonstrate 
again the advantage of keeping track of the exact switch- 
ing time. If the switching time is not known (up to 
the pulse duration r, which we take to be much longer 
than 1/E), one must use the matrices Ulg{T)Uns{T) and 
U^{0,t)Us{0,t) as the measurement matrices. The two 
eigenvalues of each one of the above matrices are equal. 
Therefore, the fidelity vanishes. This result is in agree- 
ment with the generally accepted rule that this type of 
measurement cannot be used when the qubit is biased at 
the so-called degeneracy point during the measurement. 
This example shows clearly that additional information 
can be obtained by keeping track of the exact switching 
time, provided this switching time can be determined 
with accuracy smaller than 1/E. 
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FIG. 5: The S curves (i.e. switching probabihtics as functions 
of readout bias current with fixed pulse duration) for three 
different types of detectors: a strongly coupled detector [(a); 
see Sec. IV. B], a weakly coupled fast-decohering detector [(b); 
see Sec. IV. C] and a weakly coupled slow-decohering detec- 
tor [(c); see Sec.V]. The dashed lines are the S curves for the 
qubit's clockwise- and counterclockwise-current states, i.e. the 
ones that arc obtained when the qubit is biased far from the 
degeneracy point. In the example plotted in this figure, the 
probabilities for these S curves (i.e. the dashed lines) are given 
by P = 1 - e-T^/«* with -ya = e^'"'"^' and 7ijt = e^"^ (a; is 
a dimensionless parameter that represents the bias current 
used in the readout pulse, and Jl/b, represent the switch- 
ing rates for the two different persistent-current states). The 
solid lines arc the S curves for the qubit's energy cigenstates 
when these states have the form |0) = Vol\L) + VO^\R) 
and |1) = \/0.3 \L) — \/0.7 \R). In (a) the switching proba- 
bilities for the states |L) and \R) are averaged (in a weighted 
manner) in order to obtain the switching probabilities for the 
states |0) and |I). This reasoning is the one that would be 
consistent with the picture of a strong measurement. In (b) 
the switching rates are averaged, and the switching probar 
bilities are calculated from the averaged switching rates (see 
Sec. IV.C). In (c) the parameter x is averaged, and from the 
two weighted averages of x the switching rates and probabil- 
ities are calculated (see Sec. V). Here we are treating x as 
the physical observable that is modified by the qubit at the 
Hamiltonian level. 
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FIG. 6: The maximum measurement fidelity (i.e. at the op- 
timal value of x) as a function of the angle /3 (between the 
qubit Hamiltonian's axis and the probed operator axis) for 
the three different types of detectors whose S curves are plot- 
ted in Fig. 5. The fidelity is derived by finding the maximum 
vertical separation between the two S-curves for a given value 
of /3. In (a) the fidelity is given by cos/3, as expected for the 
textbook-type projective measurement. In (b) the fidelity is 
given by Eq. (I51|l . In (c) the fidelity is calculated numeri- 
cally. For further demonstration that the coherent detector's 
fidelity can remain close to 100% even for /3 7^ 0, we plot 
in Fig. 6c (dashed line) the fidelity when the coefficient 5 in 
the switching-rate formulae of Fig. 5 is replaced by the coef- 
ficient 10 (i.e. in going from the solid to the dashed line, the 
sensitivity of the detector is increased). 



